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Abstract The universe with adiabatic matter creation 
is considered. It is thought that the negative pressure 
caused by matter creation can play the role of a dark 
energy component, and drive the accelerating expan¬ 
sion of the universe. Using the Type la supernovae (SNe 
la) data, the observational Hubble parameter data, the 
Cosmic Microwave Background (CMB) data and the 
Baryonic Acoustic Oscillation (BAO) data, we make 
constraints on the cosmological parameters, assuming a 
spatially flat universe. Our results show that the model 
with matter creation is consistent with the SNe la data, 
while the joint constraints of all these observational 
data disfavor this model. If the cosmological constant is 
taken into account, a traditional model without matter 
creation is favored by the joint observations. 
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1 Introduction 

A great encouraging development in modern cosmology 
is the discovery of the acceierating expansion of the uni¬ 
verse through observations of distant Type la Super¬ 
novae m- The anisotropy of the cosmic microwave 
background (CMB) resuits from bailoon and ground 
experiments [31[11[30[71[S] and recent WMAP m obser¬ 
vation confirmed the resuit from SNe la and favored a 
spatiaily fiat universe. It is weil known that ali known 
types of matter with positive pressure generate attrac¬ 
tive forces and deceierate the expansion of the universe. 
The discovery from SNe la and CMB indicate the exis¬ 
tence of a new component with negative pressure, which 
is now generaiiy cailed dark energy. 

The simpiest one of dark energy is a cosmoiogicai 
constant [imniiii]. An expianation of the cosmoiogicai 
constant is the vacuum energy, however, it is 120 orders 
of magnitude smaiier than the naive expectation from 
quantum fieid theory. Bothering physicists much, other 
types of dark energy are proposed, such as quintessence 
[niiiiiiii] , which is described in terms of a cosmic scaiar 
fieid (j)-, or other modified cosmoiogicai models are dis¬ 
cussed, such as the Cardassian expansion model which 
investigates the acceleration of the universe by a mod¬ 
ification to the Friedmann equations [HUH], the brane 
world model which explain the acceleration through the 
fact that general relativity is formulated in 5 dimensions 
instead of the usual 4 [IE[Ill2i| and so on. 

All of these dark energy cosmological models are 
based on the Big Bang cosmology. A model with adi¬ 
abatic matter creation was proposed firstly in order to 
interpret the cosmological entropy and solve the big- 
bang singularity problem m- The basic idea is to mod¬ 
ify the usual energy conservation law in open system 
in the framework of cosmology, which adds a balance 
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equation for the number density of the created parti¬ 
cles to the dynamic equations of the universe. Neverthe¬ 
less, after the discovery of the accelerating expansion of 
the universe, this model was reconsidered to explain it 
and got some unexpected results. The matter creation 
pressure Pc, which is negative as pointed out several 
decades ago by Zel’dovich [^, might play the role of 
a dark energy component and lead to the accelerating 
expansion of the universe. Lima & Alcaniz tested the 
model without cosmological constant through the look- 
back time-redshift relation, luminosity distance-redshift 
relation, angular size-redshift relation and the galaxy 
number counts-redshift relation Em]. It was shown 
that this model was consistent with the observational 
accelerating expansion of the universe, and could also 
alleviate the conflict between the age of the universe 
and the age of the oldest globular clusters. Zimdahl et 
al. employed the SNe la data to test the matter cre¬ 
ation scenario and also got the result of accelerating 
expansion ESj. Freaza et ah, however, based on the ob¬ 
servational SNe la data and the simulated Supernova 
Acceleration Probe (SNAP) data, showed the matter 
creation mechanism was not favored to explain the cos¬ 
mic acceleration [55]. 

In this paper, we use recently 186 SNe la sample 
[27] . combined with the observational H{z) data from 
the differential age measurements of galaxies [28], the 
CMB and BAO data [551133] . to test the cosmological 
model with matter creation and make constraints of 
the parameters. As a comparison, the model with both 
matter creation and cosmological constant is also ex¬ 
amined. This paper is organized as follows: we present 
the basic cosmological equations of the universe with 
adiabatic matter creation in Sec. [5] In Sec. [3] we give 
a brief introduction to the observational data, and the 
results and discussion are given in Sec. [4] 


2 The Cosmological Basic Equations with 
Adiabatic Matter Creation 


The Robertson-Walker(RW) metric describing the space- 
time of the universe is 

ds^ = -dt^ -I- c?{t)[zr-—r^ + -I- sin^ 6<d(/)^], (1) 

i. — kt 


where r, 0, and (j) are dimensionless comoving coordi¬ 
nates, k = Q, ±1 represent the curvature of the spatial 
section and a(t) is the scale factor. Using the Einstein 
field equation, we can acquire the equations to describe 
the dynamic behavior of the universe, namely the Fried¬ 
mann equations 


0? SttG a k 


d AttG , M A 

-=- ^[p + Hp + Pc)] + (3) 

where p = pM + Pr is the energy density(matter and 
radiation), p and A are the thermal pressure and cosmo¬ 
logical constant respectively. Pc is the matter creation 
pressure that takes the following form [34] 


Pc = - 


P + P 
SnH 




(4) 


where H = a/a is the Hubble parameter and is the 
matter creation rate. In models with adiabatic creation, 
the balance equation for the particle number density n 

is [551155] 


n „ d lb 

-h 3— = —. 

nan 

We take the form of matter creation rate as [55] 


(5) 


tb = SPnH, ( 6 ) 

where the parameter (3 is defined on the interval [0,1] 
and assumed to be constant. Matter and radiation cor¬ 
respond to /3 m and /3/j respectively. Together with the 
equation of state (EOS) 


P = wp, (7) 

the equation system becomes complete. 

From Eqs.(I5])and(15]), we can get 

d(pa^) 

Combining Eqs.(l4]),([7])and(l8]), it is easy to find 

w = 0, Pm = for matter 

w = 1/3, ph = for radiation 

ui = — I, PA = PAo for cosmological 

constant 


( 8 ) 


.(9) 


Using Eq.® and noting that 1 + z = ag/a, we can 
rewrite Eq. ([5]) in the form of the Hubble parameter 

H^=H^E^z), (10) 


where Hq is the Hubble constant, and 
E‘^{z) = f2M{l + + Hiiil + 2)4(1-^^^) 

+ (H) 


represents the expansion rate, in which [2m,^Rt^a 
and f2k are the matter density, radiation density, cos¬ 
mological constant and spatial curvature parameters at 
present. 

The deceleration parameter q is defined as 


aa 


( 12 ) 


Using Eqs.(15]),([5]) and (O, the deceleration parameter q 
reduces to 

+ (1 - 2/3fl)CH(l + z)4(i-/3«) _ 


( 2 ) 


(13) 
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In the following discussion, we consider a spatially flat 
universe [9], and neglect the radiation term for its ex¬ 
tremely small value today, that is, Qk = = 0, so 

f^M + = 1- 

If there is no cosmological constant, that is, = 
0, I7 m = 1, Eq. (ED can be rewritten as 
£;2(^) = (l + z)3(i-/5), (14) 

where [3 is /3m, and q{z) can be simplified as 

9 ( 2 ) = —^—> (15) 

which shows that q is independent of redshift. From 
Eg. lfTSl) we know that, if the universe is accelerating 
expanding, that is, qo < 0, P needs to be greater than 
1/3. 


3 The Observational Data 


SNe la are thought to be standard candles and can 
be used as distance probes. The theoretical prediction 
for luminosity distance of an astronomical object in a 
spatially flat universe is 


dL = 


c(l -I- z) 
Ho 


dz 

wr 


(16) 


where E{z) is defined in Eg. lfTUl) . The distance modulus 
is 

dz 


m = 51g^+ 25 = 42.38+ 51g[i^ [ 
Mpc h Jq 




r](17) 


where h = Ho/{10Qkm/s/Mpc). The parameter for 
SNe la is 


xlv 


= E 


[pthjzt; h, J?m,/ 3) - MSAf(+)]^ 

(TSN{zi) 


(18) 


where psNizi) is the observed distance modulus of the 
SN with redshift Zi and usN{zi) is the observational 
error. 

We also include the observational Hubble parameter 
H{z) in this work. The Hubble parameter depends on 
the differential age of the universe in this form 


H{z) 


1 dz 
1 + z dt ’ 


(19) 


which provides a direct measurement for H (z) through 
a determination of dz/dt. Using the differential ages of 
passively evolving galaxies determined from the Gem¬ 
ini Deep Deep Survey (GDDS), Simon et al. determined 
a set of observational H{z) data in the redshift range 
0 ~ 1.8 [28] . These observational H{z) data have been 
used to constrain the dark energy potential and its red¬ 
shift dependence by Simon et al. E). Various cosmo¬ 
logical models were tested using this data set in recent 
several years. Yi & Zhang used them to analyze the 


holographic dark energy models and showed that the 
fitting results did not conflict with some other indepen¬ 
dent cosmological tests [29]. Samushia & Ratra used 
the data set to constrain the parameters of HCDM, 
XCDM and i/CDM models. The constraints are consis¬ 
tent with those derived from SNe la [HD] ■ Wei & Zhang 
tested a series of other cosmological models with inter¬ 
action between dark matter and dark energy using this 
data set. m- In one word, the observational H{z) data 
are demonstrated to be an effective complementarity to 
other cosmological probes. 

For the Hubble parameter H{z) data, we have 

2 [HthiZi', h, T P) ~ 3?oi,(+)] /'on^ 

Xh — / ^ -7—N-• \Z'J) 


The model-independent shift parameter R, which 
can be derived from CMB data, is defined as 



( 21 ) 


with Zr = 1089 the redshift of recombination. From 
the three-year result of WMAP, Wang & Mukherjee 
estimated R = 1.70 ± 0.03 [32]. 

Using a large spectroscopic sample of luminous red 
galaxies from Sloan Digital Sky Survey (SDSS), Eisen- 
stein et al. successfully found the acoustic peaks in the 
bayonic matter anisotropy power spectrum, described 
by the model-independent A parameter 


A = \/ 17m 


1 

dz ■ 

[ziUl/2(^,) j 

0 E(z)_ 


12/3 


( 22 ) 


with zi = 0.35 the redshift at which the acoustic scale 
has been measured [33]. Eisenstein et al. suggested the 
value of A parameter as A = 0.469 ± 0.017. 

The model parameters h, Qm and /3 can be deter¬ 
mined through the minimization method. The com¬ 
bined can be written as 


2 2 I 2 I 

X — XsN + Xh + 


(R- 1.70)2 

0(032 


{A - 0.469)2 

0.0172 


(23) 


By minimizing x^ we can get the best fitting values of 
the parameters. 


4 Results and Discussion 

Using the observations of SNe la, Hubble parameter 
i7(z), CMB and BAG, we test the cosmological model 
with adiabatic matter creation, assuming a spatially flat 
universe. We show in Eig. |T]the confidence regions of 
parameters h and /3 of this model, for different sets of 
observational data. The best fitting results of parame¬ 
ters are listed in Table [TJ We can see from this figure 
that, for the SNe la data, /3 is greater that 1/3 at 3cr 
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Fig. 1 Confidence regions of 1, 2 and 3 a (from inside to outside) 
for the two parameters h and f3, for different observational data 
set as labelled in the figure. The horizon lines are of = 1/3 and 
0 . 


Fig. 2 Confidence regions for the two parameters Hm and /9,for 
the cosmological model with both the matter creation and the 
cosmological constant. Parameter h has been marginalized. The 
lines are coded as Fig. [T] 


Table 1 Best fitting values of h and f} of the matter creation 
model without the cosmological constant. 

Test h f3 

SNe la 0.64 0.54 239.97/184 

H(z) 0.66 0.33 8.50/7 

SNe Ia+H(z) 0.62 0.42 269.40/193 

SNe Ia+H(z)+CMB+BAO 0.56 -0.04 753.71/195 


Table 2 Best fitting values of Hm and /3 of the matter creation 
model with the cosmological constant {h marginalized). 

Test 17 m /3 x^/d.o.f 

SNe la 0.10 -0.66 237.64/184 

H(z) 0.91 0.31 13.94/7 

SNe Ia+H(z) 0.19 -0.26 250.87/193 

SNe Ia+H(z)+CMB+BAO 0.30 -0.02 256.42/195 


level, which means that the universe is accelerating ex¬ 
panding. The results of ll{z) data are a little different 
from which of SNe la. The best htting value of /3 is 
0.33 for li{z) data. Because of the smallness of statis¬ 
tical sample of H{z) data, the error bars of parameters 
are relatively large and the H{z) data fail to constrain 
the parameters strictly enough. H{z) data cannot con¬ 
strain the model parameters very well. The joint con¬ 
straints of all these data (SNe la -I- H(z) -I- CMB -I- 
BAO) are plotted in the same figure. We find f3 < 0 
at 3cr level, with the best fitting value jS = —0.04. A 
negative value of [3 means matter vanishing instead of 
matter creation, which is forbidden thermodynamically 
[M] . We draw the conclusion that this model is consis¬ 
tent with SNe la observations, but is disfavored by the 
joint analysis from different observations. 

One can expect that the negative pressure of both 
the cosmological constant and the matter creation may 


jointly generate the acceleration of the universe. We 
also test the model with a cosmological constant in or¬ 
der to give a comparison. After marginalizing parame¬ 
ter h, we show the confidence regions of parameters f^M 
and /3 in Fig. [5] and the best fitting parameters in Table 
[2] respectively. The cosmological constant Ha = 1 —17m- 
For the SNe la and H{z) data, the constraints of model 
parameters are not very strong. The statistical uncer¬ 
tainties of model parameters are so large that we cannot 
give any convincible conclusions. While the joint con¬ 
straints of these 4 kinds of observations show that [3 
tends to be zero, which reduces to the familar ACDM 
cosmology. The best fitting result 17 m = 0.30 ± 0.02 
is also consistent with other studies of the traditional 
ACDM model [27ll^. We also notice that our result 
is very similar with [26j . However, in their treatment, 
a small SNe la sample (16 low-redshift and 38 high- 
redshift supernovae from Perlmutter et al. m) could not 
give strong constrains on the model parameters. They 
drew the conclusion by adopting a Gaussian prior of 
17 m = 0.27 ± 0.06. In our work, we give the constraints 
from joint astronomical observations directly. 

In summary, the idea that a negative matter cre¬ 
ation pressure may play the role of dark energy and 
drive the accelerating expansion of universe, seems not 
favored by the combined observations. While for the 
model with both the adiabatic matter creation and cos¬ 
mological constant, the joint constraints tend to give 
/3 ~ 0, which means no matter creation, and the model 
becomes the traditional ACDM one. 

Acknowledgements We would like to thank Yue Bin, Dr.Peng- 
Jie Zhang, Li Chen and Zong-Hong Zhu, for their friendly discus¬ 
sion. This work was supported by the National Science Founda¬ 
tion of China (Grants No. 10473002 and 10273003), the Scientific 















5 


Research Foundation for the Returned Overseas Chinese Schol¬ 
ars, State Education Ministry and the Research Foundation for 
Undergraduate of Beijing Normal University. 


References 

1. Riess, A.G., Filippenko, A.V., Challis, P., et al.:AJ 116, 
1009(1998) 

2. Perlmutter, S., Aldering, G., Goldhaber, G., et al.:ApJ 517, 
565(1999) 

3. Miller, A.D., Galdwell, R., Devlin, M.J., et al.:ApJ 524, 
Ll(1999) 

4. de Bernardis, P., Ade, P.A.R., Bock, J.J., et al.:Nature 404, 
955(2000) 

5. Hanany, S., Ade, P., Balbi, A., et al.:ApJ 545, L5(2000) 

6. Halverson, N.W., Leitch, E.M., Pryke, G.:ApJ 568, 38(2002) 

7. Mason, B.S., Pearson, T.J., Readliead, A.G.S., et al.:ApJ 591, 
540(2003) 

8. Benoit, A., Ade, P., Amblard, A., et al.:A&A 399, L25(2003) 

9. Spergel, D.N., Verde, L., Peiris, H.V., et al.:ApJS 148, 
175(2003) 

10. Weinberg, S.:Rev. Mod. Pliys. 61, 1(1989) 

11. Carroll, S.M., Press, W.H., Turner, E.L.:ARA&:A 30, 
499(1992) 

12. Ostriker, J.P., Steinhardt, P.J.:Nature 377, 600(1999) 

13. Ratra, B., Peebles, P.J.E.iPliys. Rev. D. 37, 3406(1999) 

14. Coble, K., Dodelson, S., Frieman, J.A.:Phys. Rev. D. 55, 
1851(1997) 

15. Caldwell, R.R., Dave, R., Steinhardt, P.J.:Phys. Rev. Lett. 
80, 1582(1998) 

16. Freese, K., Lewis, M.:Phys. Lett. B. 540, 1(2002) 

17. Zhu, Z., Fujimoto, M.:ApJ 581, 1(2002) 

18. Randall, L., Sundrum, R.:Phys. Rev. Lett. 83, 3370(1999) 

19. DefTayet, C., Dvali, G., Gabadadze, G.:Phys. Rev. D. 65, 
044023(2002) 

20. Avelino, P.P., Martins, C.J.A.P.:ApJ 565, 661(2002) 

21. Prigogine, I., Geheniau, J., Gunzig, E., Nardone, P.:Gen. Rel. 
Grav. 21, 767(1989) 

22. Zel’dovich, Ya. B.ijournal of Experimental and Theoretical 
Physics 12, 307(1970) 

23. Lima, J.A.S., Alcaniz, J.S.:A&A 348, 1(1999) 

24. Alcaniz, J.S., Lima, J.A.S.:A&:A 349, 729(1999) 

25. Zimdahl, W., Schwarz, D.J., Balakin, A.B., Pavon, D.:Phys. 
Rev. D. 64, 063501(2001) 

26. Freaza, M.P., de Souza, R.S., Waga, I.:Phys. Rev. D. 66, 
103502(2002) 

27. Riess, A.G., Strolger, L., Tonry, J., et al.:ApJ 607, 665(2004) 

28. Simon, J., Verde, L., Jimenez, R.:Phys. Rev. D. 71, 
123001(2005) 

29. Yi, Z.L., Zhang, T.J.:Mod. Phys. Lett. A. 22, 41(2007) 

30. Samushia, L., Ratra, B.: ApJ 650, L5(2006) 

31. Wei, H., Zhang, S.N.:Phys. Lett. B. 644, 7(2007) 

32. Wang, Y., Mukherjee, P.:ApJ 650, 1(2006) 

33. Eisenstein, D.J., Zehavi, I., Hogg, D.W., et al.:ApJ 633, 
560(2005) 

34. Calvao, M.O., Lima, J.A.S., Waga, L:Phys. Lett. A. 162, 
233(1992) 

35. Lima, J.A.S., Germano, A.S.M.iPhys. Lett. A. 170, 
373(1992) 




